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We develop an action principle for those models arising from isotropic loop quantum cosmology, 
and show that there is a natural conserved quantity Q for the discrete difference equation arising 
from the Hamiltonian constraint. This quantity Q relates the semi-classical limit of the wavefunction 
at large values of the spatial volume, but opposite triad orientations. Moreover, there is a similar 
quantity for generic difference equations of one parameter arising from a self-adjoint operator. 
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I. INTRODUCTION 

Loop quantum gravity (LQG) has been a successful 
research program, in the sense that it is a rigorously de- 
rived picture of the behavior of the gravitational field 
in the quantum regime. However, it is difficult to ap- 
ply as is to physical problems - such as how quantum 
gravity effects in the realm of classical singularities af- 
fect cosmological models in general - because of its com- 
plexity. This has led to the development of loop quan- 
tum cosmology (LQC), a symmetry-reduced version of 
LQG. The Hamiltonian constraint equations of LQC are 
fundamentally discrete, due to the corresponding nature 
of quantum geometry coming from LQG, and thus are 
difference equations, rather than differential equations. 
The study of such difference equations has allowed the 
use of various analytic solution techniques, e.g. gener- 
ating functions but these quickly become difficult to 
work with, leaving only numerical methods [2|]. Thus, 
it is fruitful to look for other analytic methods to ob- 
tain exact information about solutions to the quantum 
constraint equation; we do so by considering the discrete 
analogue of classical mechanics. 

Discrete Lagrangian mechanics has been developed in- 
dependently several times, in many contexts; for a list of 
previous work, see Marsden and West |3J • Specifically for 
LQC, these ideas were developed by Shojai and Shojai Q, 
but oriented more towards finding approximate solutions 
to difference equations. There is also the program of 
consistent discretizations [f|, aiming to solve discretized 
versions of constrained mechanics, so that the constraints 
are preserved under evolution of the system in the model 
parameter n (for an isotropic model). There, the em- 
phasis is on the algebraic consistency of the discretized 
mechanics. In this work, we discuss a discrete version of 
the familiar Lagrangian mechanics, and show that for the 
isotropic models currently used in LQC, there is a con- 
served quantity Q - in the sense that this function Q n of 
the single (discrete) parameter n arising in the quantum 
Hamiltonian constraint equation is constant, regardless 
of what n is evaluated at. 



This paper is organized as follows. We start with an 
action principle for a functional C n of a one-parameter 
sequence s n , and derive the discrete version of the Eulcr- 
Lagrange equation which extremizes the action. Then, 
assuming the action is invariant under an infinitesimal 
transformation of s n , we find a quantity Q that re- 
mains constant under the action of the difference equa- 
tion for s n . As we will show, the generic Lagrangian 
for the isotropic models will be phase invariant under 
s n — » s' n = s n exp(ia); this generates a discrete version of 
the familiar conserved quantity for a phase invariant field 
4>(t). After this groundwork is laid, we turn to specific ex- 
amples in isotropic LQC, and find the corresponding con- 
served quantities Q. We also consider general self-adjoint 
difference operators (or difference equations of the same 
form) and show they similarly have a conserved quantity. 



II. AN ACTION PRINCIPLE FOR 
DIFFERENCE EQUATIONS 

For the rest of this work, we define the difference op- 
erator A acting on a sequence s n as 



As n 3 n Sji—X, 



(1) 



One particularly useful identity is the discrete version of 
integration by parts, namely the equivalence 



M M 

F n AG n = 



[A(F n+1 G n ) - G„AF n+1 }. 

(2) 

Here the summation acts as the discrete analogue of a 
integral for functions of continuous variables. Much like 
the integral of a total derivative becomes a boundary 
term, the sum of a total difference gives a similar result: 



M 



A(F n+ iG n ) — Fm+iGm — F-m+iG- 



M ■ 



(3) 



n=-M+l 
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The "unbalanced" limits in these series come from the 
definition of A. These limits have also been chosen to be 
symmetric around n = 0, but this is not a requirement; 
all of the results obtained in this work would be similar 
for a different choice of limits. 
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Now we briefly derive the Euler-Lagrange equation for 
a one-parameter sequence, both to keep the discussion 
self-contained, and to develop the notation used here. 
We start with the action 



M 

E 

n=— Af+1 



n= -M+l 



The format of C n is to match that of "standard" continu- 
ous actions, i.e. with a coordinate variable q and its time 
derivative q. We will assume that the boundary values 
sm and are fixed data, and determine the rest of 
the sequence. First, we seek the general Euler-Lagrange 
equation for the sequence s n that extremizes the value of 
S. In this vein, we consider an infinitesimal variation of 
the sequence s n , namely, 



Sn -> S,,. 



for | tt. | < M, and find the sequence such that 

M 



dS_ 

da 

This gives 



E 



dC{n;s' n , As' n ) 



n= -M+l 



da 



8C n 



0. 



a=0 



da ds n ™ d(As n ) 



(At7„). 



(4) 



When placed inside the summation, the last term in the 
equation (@| can be integrated by parts as in J2J), i.e. 



M 



E 

n=-M+l 
M 



dC(n; s n , As n ) 



(Affc, 



= E a 

n=— Af+1 
M 

- E a 

n=-M+l 



d(As n ) 
d£(n+ 1; s n+ i, As n+1 ) 



d(As n+1 ) 
dC{n+ l;s„+i, As n+ i 



->)n 



d(As n +i 

From equation ([3]), the total difference in ([5]) becomes 
~d£(n + l;s n+1 ,As n+ x) 



(5) 



M 

E A 

n=-M+l 



d(As n+ i) 
dC(n;s n , As n ) 



-'In 



(G) 



d(As n ) 
dC(n; s n , As n ) 



Vm 



n=M+l 



d(As 



V-M = 0. 



n=-M+l 



where the last equality results from the assumption that 
the boundary values are fixed, so t\m = V-m — 0. Thus, 



clS_ 

da 



M-l 

E 



dC{n;s n , As n ) 
ds n 



n=-M+l 

a / dC(n + l;s n+ i,As n+ i) 
^ d(As n+1 ) 



The remaining values r\ n are arbitrary, so this gives us 
the Euler-Lagrange equation^] 



dC(n;s n ,As 7 
ds n 



dC(n+ l;s„ +1 , As n+1 ) 



d(As n+1 ) 



= 0, (7) 



for \n\ < M - 1. 

As an example that will be relevant later, suppose we 
consider the action 



S 



M 

E 

= -M+l 



(f n \As n \ 



■9n\Sn\ 



(8) 



where s n and its complex conjugate s n are considered 
independent - thus there are two equations of motion, 
one for each sequence - and the coefficient functions f n 
and g n are real. Then the Euler-Lagrange equations are 

g n s n - A(/„ + iAs n+ i) = 

(and conjugate), or 

fn+lSn+l — (fn+1 + fn) s n + fn s n-l = 9nS n - (9) 

As we will see below, this is the same form as the equation 
for the eigensequences in the LQC isotropic model. For 
general equations of this form, there are two solutions, 
since the initial data can be independently selected at 
two values of n. 



III. CONSERVED QUANTITIES 

Suppose that the Lagrangian C(n; s n , As n ) is invariant 
under a transformation of the sequence s„, i.e. 

&ti ^ s n -\- a^ ni 

but now including variations of the boundary values Sm 
and s_m! this will give us a conserved charge Q n , as a 
special case of the discrete Noether's theorem. Since this 
is essentially what we have above in equation ([5]), then 
using the Euler-Lagrange equations, 



M 

E a 

n=-M+l 



dC(n + 1; s n+1 , As n+ i; 
9(As n+ i) 



0. (10) 



1 We note here that another form of this equation has been derived 
in the literature. If the original discrete Lagrangian is of the form 
C(n; s n , s n _i) - that is, all appearances of s n _i in C n are not 
necessarily part of a difference As n , then the Euler-Lagrange 
equations are of the form 

dC(n; l) 8C{n + 1; s n+ i,s n ) _ ^ 

ds n ds„ 
for \n\ < M — 1. Note that the position of the variables 
is important in the above equation, since L>ri 

(u, v) ^ C„(v, if). 

This is the form usually used in the literature, such as by Bahr, 
Gambini and Pullin and Shojai and Shojai although the 
latter allows the Lagrangian to be a function of s n -2 as well as 
s n ,s n —i. If £ is a function of s n and As n only, however, this is 
equivalent to the equation (0. 
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One way to interpret this is that a charge Q n , given by 
n _ dC(n + 1; s n+ i, As n+1 ) 

W — a/ a \ Sin 

so that AQ n = 0. Using the Euler-Lagrange equa- 
tions, we can show that AQ n = for all n, not just the 
sum over the entire range. Therefore we have Q n = Q = 
constant for all |n| < M. 

We turn again now to the example action ([8]). which is 
invariant under s n — > s' n = s„exp(ia). Therefore, in the 
infinitesimal limit, £ n = is n . Since the sequences s n and 
s„ are considered independent, the current is 

dC n +i 



_ dC n+ i 

^ ~ A \ V 



<9(As n+ i) d{As n+1 ) 
ifn+i [s n {As n+1 ) - s n (As n+ i)] 



(11) 



Here, we make two comments. The first is that Q de- 
pends only on the "kinetic" coefficient function /„, and 
not the "potential" function g n . Thus if two cosmological 
models have constraint equations where only the func- 
tion g n differs between the two, they will have exactly 
the same form of conserved quantity Q. This occurs for 
Friedmann-Robertson- Walker models, where the value of 
the curvature constant k changes only the functional form 
of g n . Finally, we note this is a discrete version of the 
conserved current for a complex scalar field. In particu- 
lar, in the large parameter regime |n| 3> 1, we can use 
a Taylor series expansion of s rl +i and s n +i, so that the 
sequence is approximated by a continuous function s(n). 
This limit gives 



ifo(n) 



s(n) 



ds(n) 
dn 



s(n) 



ds(n) 



dn 



(12) 



where fo(n) is the lowest order term in the series expan- 
sion of fn+i- In the context of a fundamentally discrete 
quantum cosmology, this Q would be the semi-classical 
limit of this conserved quantity. 

Before we turn to specific cases in isotropic LQC, we 
point out that a similar analysis is possible for differ- 
ence equations with multiple parameters n*, rather than 
the single parameter n used here. However, this exten- 
sion would require the difference equation to have con- 
stant lattice spacing. Such a constraint was developed 
for Bianchi I (6|, although there is some debate whether 
physical considerations rule out this model in lieu of a 
difference equation with lattice refinement, i.e. the step 
sizes depend on the parameters n, themselves Q • In any 
case, the derivations presented above cannot be extended 
to such lattice refined models. 



IV. QUANTIZATION OF LQC ISOTROPIC 
MODEL 



For the first use of the general methods developed 
above, we consider the flat isotropic LQC model of 



Ashtekar, Pawlowski and Singh (APS) [8]. In this model, 
they consider a scalar field on a Friedman-Robertson- 
Walker space-time, and use this matter field as the in- 
ternal clock for evolution. The sequences themselves are 
parametrized by a value proportional to the eigen- 

value of the volume operator; the sign of v gives the 
orientation of the spatial triad. Using the notation of 
APS, a sequence of values {^(v)} that solves the overall 
Hamiltonian constraint can be decomposed into an inte- 
gral over eigensequences e u (v) with eigenvalues u>, solving 
the equation 

C+iv^iv + 4) + C»e w ( W ) + C-(v)e u (v - 4) 
= oj 2 B(v)eUv), (13) 



where 



C+(v) = -^— \v + 2\\\v + 1|-|« + 3||, (14a) 



C-(v)=C + (v-4), 
C°(v) = -C+(v)-C-(v). 



(14b) 
(14c) 



and 



B(v) 



2 K ''' 



lu + llVa.^.iji/s 



with the constant K = 2v / 3/(3\/3v / 3). The solutions 
e w (t>) of the difference equation (|13l) have support only 
on values v = 4n + e, where e € [0,4). With this in mind, 
the relation (ITBl can be put into the Euler-Lagrange form 
© if we let 

v = 4n + e, 

for e G [0,4) and 

s n = e u (v), f n = C~(v), g n = uj 2 B(v). 

Thus, these eigensequences can be derived from a La- 
grangian of the form 

C v = C~ (v) | Ae w (v) | 2 + lo 2 B{v) \e u (v) | 2 , 

which has a conserved quantity 

Q = iC + ' [v)[e u {v)e u [v + 4) - e u {v + 4)e u (v)]. (15) 

In particular, note that C + (0) ^ 0, for those wave func- 
tions having support at the classical singularity v = 0; 
this means that the conserved quantity Q is not required 
to be zero for those solutions. This allows the use of this 
quantity to relate the asymptotic linear combinations of 
the wave function at large values of the volume, which 
we examine next. 

In the semi-classical limit of large spatial volumes 
\v\ » 1, we expect to recover the usual Wheeler-deWitt 
(WdW) theory for this space-time. Specifically, tak- 
ing the Taylor series of the difference equation (fT3| at 
large v gives a differential equation which matches the 
WdW equation at lowest order. These WdW solutions 
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are worked out by APS, and the eigenfunctions are given 
by 



e ±\k\{v) 



1 



,±»|fc|ln|u| 



(16) 



where ui = | k \ yYhiG. When \v\ ^> 1, the two indepen- 
dent solutions e u to the isotropic LQC difference equa- 
tion, for a given uj, in the semi-classical limit should each 
match a linear combination of the solutions e±|jw to the 
corresponding WdW equation: 

e u (v) -> Ae\ k \(v) + Be_\ k \(v) v > 1, 
e u (v) -> Ce| fe |(u) + De_\ k \(v) v -C -1, 

for constant coefficients {A, _B, C, I?}. In the large vol- 
ume limit, the conserved current becomes 



Q 



SiirKG 
16 1 



e w («) 



dv 



dp 

e u (v) — 
dv 



(17) 



This asymptotic limit of Q is precisely the same as the 
analogous conserved quantity that is obtained in WdW 
theory, using Noether's theorem for the appropriate con- 
tinuous action. 

From (fTB)) we have 

deu ikA ikB 

ir^^^-^iw- (18) 

In the large positive volume limit - i.e. large volume for 
a positively oriented triad - (fTT)) gives 

Q™ = ^f^(L4| 2 -|3| 2 ), 

while in the large negative volume limit (large volume 
but negative orientation of the triad), 



Qn- 



Z-rtkKG 
16 



{\C\ 2 -\Df). 



Thus, the conservation of the current (fTT)) gives 
\A\ 2 - \B\ 2 = \C\ 2 - \D\ 2 



(19) 



as shown numericalljB for their quantization Q. Note 
that the coefficient function g n does not enter at all into 
the form of the conserved quantity Q. Indeed, since go- 
ing from the k = isotropic model to either k = ±1 
only alters this function g n |10| , the discussion here goes 
through without change. 



2 In APS, the authors choose to look only at those quantum solu- 
tions of the discrete LQC Hamiltonian constraint operator that 
are symmetric under the parity operator, to match with the 
choice of WdW eigenfunctions (1 1 6 D . However, these sequences 
are constructed out of linear combinations of the eigensequences 
e u , which are not symmetric under the exchange v — > —v. Thus, 
the conservation law J 19H is non-trivial. 



V. SELF-ADJOINT HAMILTONIAN 
CONSTRAINT OPERATORS 

The last section dealt with a particular choice of factor- 
ordering for the Hamiltonian constraint of isotropic LQC, 
so it is interesting to see how generic these results are. 
The particular property we focus on here is the self- 
adjoint nature of the constraint operator; this is neces- 
sary in APS in order to find Dirac observables or physical 
states, chosen with a Hilbert space structure. Thus, in 
the following, we look for a conserved quantity Q for a 
self-adjoint (gravitational) Hamiltonian constraint oper- 
ator H g acting on states ip. 

For the isotropic model, we have a basis of states \/i) 
which are eigenstates of the volume operator, where \x € 
M. There is a natural inner product for these states, 
namely, 

(n\v) = 5^ v . 

Note this is a Kronecker delta, not a Dirac delta function. 
Suppose we have an operator A acting on basis states \fj) 
such that 

= A+((x)\» + S) + A°(ji)\n) + A-(ji)\n - S). 

where 8 is a constant value indicating the step size for 
the lattice used in the difference equation; the value of 5 
results from the holonomy operator acting on the state 
and is related to a physical choice, such as the mini- 
mum length of the model. Then, one can show that A is 
self-adjoint, i.e. (Afi\u) — (n\Av) only if 

A+( f i) = A-(» + 8). 

If we define a physical wavefunction as ip^ = (ip\fj), then 
the requirement that the gravitational Hamiltonian con- 
straint operator H g is self-adjoint means that the con- 
straint 

Hipp = (H g + H ma , tt )iptx = 

gives rise to a difference equation of the form 

H+ifi^+sicP) + H°(m%(<£) + H-{fjL)^ 5 {4>) 
= -H matt {n)ip„(<P) (20) 

with the gravitational operator H g coefficient functions 
H°([x) and H~(fx) = H + (/.i — 8). Here, we included the 
possibility of one or more matter fields <p a matter Hamil- 
tonian operator H matt - 

For isotropic matter, we have generically that H mat t 
acts only on the matter fields <fi, giving at most a pre- 
factor B(fi) related to the volume of the space-time with 
triad eigenvalue /z. If we use the viewpoint of APS, we 
use the matter field 4> a s an internal clock in the model, 
so "evolution" means changing <p. By finding solutions 
to the difference equation 

H+i^+sicP) + H°{n)M<f>) + H-(fx)^ s (<f>) 
^ 2 B(/i)^(0). (21) 
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one can build up an arbitrary solution of the original con- 
straint (|20p out of linear combinations of eigensequence 
solutions to (|2"Tj) . To match up with the previous dis- 
cussion, one can use a new parameter n = ji/8 to get 
a lattice with unit spacing, then convert back into the 
physical parameter /x. Thus, this equation is derivable 
from an action principle based on the Lagrangian 



where the function B(fi) relates to the dependence of 
the matter operator on the metric components. For the 
constraint equation ()22j) . the step size 6 = 4/io; the pa- 
rameter /x of this function can be rescaled as n = /i/4/io 
to give unit lattice spacing. We define the new function 

Sn{(t>) = [V{in+l)vo ~ V(4n-l)fi ]'>p4nfi (<f>), (24) 

so that the constraint equation (|2"2"T) is now 



where 

F(/x) = u 2 B{^) - H°(ti - H+(p) - H~(ji) 
and conserved quantity 

Q = tH + (/i)(^^+<5 - ^M+<?VV) 

As commented above, the "potential" function is 
not relevant in finding the conserved quantity Q. There- 
fore, when the constraint operator H g is self-adjoint, 
there is a conserved quantity Q associated with the 
eigensequences of the operator, regardless of factor or- 
dering issues, although specific orderings may lead to dif- 
fering coefficient functions H + (/j,). 



VI. NON-SELF-ADJOINT CONSTRAINT 
EQUATIONS IN LQC 

We have seen that a generic self-adjoint constraint 
equation is derivable from a discrete action, with a 
resulting conserved charge from the symmetry s n — > 
s' n = s n exp(ia). However, a conserved quantity exists 
even when the original Hamiltonian operator is not self- 
adjoint, but can be written in the form (|21j) . For ex- 
ample, with the earlier quantization of the flat isotropic 
model UM, we have that 



Sn+l(^) - 2s„(0) + S„_i(0) = 



uj 2 B(n)s n ((/)) 



which is of the same form as the discrete Euler-Lagrange 
equation © when one chooses f n = 1. Note that (f2"5)l 
and (f2~4")l together requires that s — 0, so there is no 
issue with evaluating the right-hand side of this derived 
difference equation at n = 0. This difference equation 
gives a conserved quantity, which when written in terms 
of the original parameter /i is 

x [VvO)^V+4«> (0) _ VV+4 W) (0)^(0)] ( 25 ) 

The conservation of this function can be checked directly 
using the original difference equation (|22|) . 

A crucial difference between this conserved quantity 
and the charge (TT5|) obtained from the APS quantization 
is what happens to the prefactor of Q at the classical 
singularity. As stated above, for the APS model (fT"3")) . at 
this singularity, C + (0) ^ 0, so Q can take any real value; 
this allows the charge for wave functions passing through 
the classical singularity to be non-zero, and thus provide 
a relation between semi-classical limits of the wave func- 
tion far away from the v = point. On the other hand, 
for the earlier quantization (|2"2")l . we have the difference 



in volumes V u 



is 



-2(V„ +M0 - V^M*) 



87rG 7 Vo^ 2 



(22) 



where fj,o is a constant step size, related to the minimum 
eigenvalue of the area operator, 7 is the Immirzi param- 
eter, lp the Planck length, and 



TIM] 
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3/2 



(23) 



are the eigenvalues of the volume operator. The precise 
form of the operator H mat t is not important, only that it 
acts on the matter fields 4> and not directly on the triad 
eigenvalues /i. If we write the <f> dependence ip^{4>) = 
ipp exp(zu;</)), this operator will have an action of the form 



5 tlmatt^y. 



Vu, 



V- L 







so that the charge Q = for any wave function passing 
through the (i — classical singularity. This results from 
the non-self-adjoint nature of the constraint and places 
a strong restriction on allowable wave functions for this 
quantization not present for the APS model. 

As mentioned previously, sq = always, so the value 
ipo is undetermined for the same reason Q = for all 
sequences f/v ~ the fact that the coefficient function of 
4>n(4>) m the constraint equation is zero when evaluated 
at fi = 0. Thus, the equation (|22[) puts a constraint on 
the other values ip^, and thus restricts the initial data for 
the wave function. The space of solutions for this equa- 
tion is two-dimensional; thus, the restriction on ip n means 
there is a unique (up to scaling) solution to the LQC con- 
straint, a situation known as dynamical initial conditions 
- the evolution equation itself picks out the wave func- 
tion without additional physical or theoretical input [T2j . 
This has the bonus that there is no requirement of choos- 
ing a specific boundary condition for the wave function, a 



6 



choice which may lead to differing physical results. How- 
ever, in the context of discrete equations, this may lead 
to radically different behavior on either side of the clas- 
sical singularity /i — (specifically, the absence of 
pre-classical solutions for a particular orientation of the 
triad), although this issue has not been explored with 
self-adjoint constraints. 



This leads to the question of whether imposing a condi- 
tion Q = for other models may lead to a similar result; 
strictly speaking, this is not the same as the original idea 
of dynamical initial conditions, since we have added a 
specific input to the model beyond the constraint equa- 
tion, but it is worthy to look into this choice as an al- 
ternative to various boundary conditions. The answer is 
negative, which we now show for the APS quantization 
(fT5| . Here, the wave functions (TTB)) are chosen as an or- 
thonormal basis of the WdW solution space, and used 
to find the initial data for the LQC eigensequences, but 
these functions are certainly not the only choice. Indeed, 



one can pick 

h,k( v ) = + e -k(v)} = -^cos(|fc| ln|«|) 

V2 V 71 " 

i 1 
h,k{v) = -^i e \k\(v) - e_| fc |(u)] = -^=sin(|fc|ln|t>|) 

For these functions, we have that Q — 0, so any linear 
combination of them also has zero charge; any eigense- 
quence solving the quantum constraint (|13[) . with initial 
data given by such a combination for \v\ ^> 1 will have 
the same zero charge. The semi-classical analysis of APS 
can carry forward from this point, writing the generic 
solution to the semi-classical differential equation using 
the eigenfunctions fi^ and /b.fc, i.e. 



As with the APS analysis, this leads to superselection, 
and one can restrict solutions to only or VE^, where ^ 
is the solution written only in terms of the eigenfunctions 
fi t h- Thus, imposing Q = docs not lead a reduction in 
the solution space. 
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